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Abstract 

We develop a methodology to ensure that the stress tensor, regardless of its number of 
independent components, can be reduced to an exactly equivalent one which has the same 
number of independent components as the surface force. It is applicable to the momentum 
balance if the shear viscosity is constant. A direct application of this method to the energy 
balance also leads to a reduction of the dissipation rate of kinetic energy. Following this 
procedure, significant saving in analysis and computation may be achieved. For turbulent 
flows, this strategy immediately implies that a given Reynolds stress model can always 
be replaced by a reduced one before putting it into computation. Furthermore, we show 
how the modeling of Reynolds stress tensor can be reduced to that of the mean turbulent 
Lamb vector alone, which is much simpler. As a first step of this alternative modeling 
development, we derive the governing equations for the Lamb vector and its square. These 
equations form a basis of new second-order closure schemes and, we believe, should be 
favorably compared to that of traditional Reynolds stress transport equation. 


* This research was supported in part by Office of Naval Research with Dr. E.P. Rood as 
technical monitor, and the National Aeronautics and Space Administration under NASA 
Contract No. NASl-19480 while the first and second authors were in residence at the 
Institute for Computer Applications in Science and Engineering (ICASE), NASA Langley 
Research Center. Hampton, VA 23681-0001. 
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1. Introduction 


Fluid dynamics needs very intensive analyses and computations, and a reduction of 
this work can be very useful. In this paper, we demonstrate that this can be achieved for the 
balance of momentum and energy, for both laminar and turbulent flows, and this reduction 
has roots in the very fundamental nature of fluid motion. Furthermore, our analysis 
suggests that the mean turbulent Lamb vector has advantages over the traditional Reynolds 
stress; its transport equation is derived as a foundation for an alternative approach in 
turbulence modeling. 

2. Reduced Stress Tensor for Newtonian Fluid 

Let S and S' be two arbitrary differentiable second-rank tensors. If their difference is 
a divergence-free tensor, there must be V • S = V • S' = F, say. We call S and S' the tensor 
potential of the vector F ; clearly, the number of potentials for F is infinite. This gives us 
a chance to construct the simplest potential, whose number of independent components 
should be the same as that of F. As a theorem, this construction is achievable through 
making the Stokes-Helmholtz (S-H) decomposition for F: 

F = Vd + V x A, V • A = 0, (1) 

where <j> and A are the scalar and vector S-H potentials of F, and have altogether three 
independent components. It is well known that decomposition (1) exists globally for any 
integrable F. Now, let I be the unit tensor and write 

S = 4>l + 1 x -A, or S tj = 4>8 tJ - e ijk A k , (2) 

then V • S is nothing but (1). This S is exactly the desired three-component tensor 
potential of F , which we refer to as the S-H tensor potential of F. Note that it consists 
of an anisotropic part and a skew-symmetric part, whereas itself is neither symmetric nor 
skew- sv mmet r i c . 

Since the S-H decomposition stands at the center of our reasoning, we make some 
remarks on it. It is well known that for a given vector F, finding its potentials <j> and A 
amounts to solving a scalar and a vector Poisson equations, of which the integral repre- 
sentation is in terms of the Green’s function of Laplace operator. However, this classical 
representation is often inconvenient in practice. If the flow is unbounded or with periodic 
boundary conditions, <j> and A can of course be easily found in the Fourier space; but the 
best approach is using the helical-wave decomposition (HWD), with the eigenvectors of 
the curl operator forming a complete orthonormal basis. In fact, HWD is nothing but a 
further sharpening of the S-H decomposition. 1 Whereas the existing theory and application 
of HWD have been confined to unbounded flow, for which the basis vectors can be easily 
obtained in Fourier space or in terms of some special functions, 1 ,2,3 it has been proven by 
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Yoshida and Giga 4 that the complete orthonormal HWD basis also exists in an arbitrary 
bounded domain. In this paper we shall not go into these details (see Wu et al . 5 ) ; rather, 
we simply assume that, in any flow domain and for any vector of our concern, the splitting 
of this vector into the longitudinal (curl-free) part and transverse (divergence-free) part 
can be obtained readily using the procedures discussed above. 

Now, consider a Newtonian fluid, of which the viscous stress tensor is proportional to 
the strain-rate tensor D = D T , where the superscript T means transpose. Let d = V • u 
and w = V x u be the dilatation and vorticity, respectively, and fl = — £l T be the anti- 
symmetric spin tensor such that Vu = D + fi. Moreover, let p and A be the first and 
second dynamic viscosities. Then, it is easily seen that there exists an identity 6-8 

D = i?I + ft-B, 

where 

B = ill — Vu T with V • B = 0 (3a, 6 ) 

is known as the surface-strain-rate tensor (Dishington 9 , whose definition differs from (3a) 
by a transpose). Therefore, from the Cauchy- Poisson constitutive equation one obtains an 
intrinsic triple decomposition of the stress tensor T: 

T = -III + 2pft - 2pB, (4a) 

where 

II = p — (A + 2p)i9 (46) 

is the isotropic part of T. Correspondingly, we have an intrinsic triple decomposition of 
the surface stress (the traction) at any surface element of unit area in the fluid or on its 
boundary: 

t = n ■ T = -Iln + r + t s , (5a) 

where —Iln is the normal stress due to compression/expansion, and 

t = pu; x n, t s = -2 pn - B (56, c) 

are the shear stress and stress due to the resistance of the viscous fluid surface to its strain, 
respectively. 

It is now evident that, as long as p. = constant, then pB is divergence-less and can well 
be dropped off from the momentum equation. What left is precisely the three-component 
S-H tensor potential 

T = -III + 2 pft, (6a) 

such that the Navier-Stokes equation with an external body force p/, 

Dti , ~ 

p— - p/ = V • T = -VII - V x (pw) (66) 
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represents a natural S-H decomposition of the total body force (inertial plus external). 

Although equation (6b), in particular its incompressible version, has been known for 
long time, 10 ' 11 the concept of reduced stress tensor (6a) is not widely appreciated, and its 
application in practice has not been fully explored. The only exception known to the authors 
is the work of Eraslan et al. 12 , who, in developing a numerical method, noticed that some 
components of T tJ played no role in their scheme and simply ignored them. A cubic mesh 
box. originally having 27 control points, was thereby simplified to having 7 points only. 
This reduced the CPU time for estimating stress by 70%, and for the overall computation, 
by 40%. This example indicates that the saving due to intentionally replacing T by T can 
be very significant. 

3. Self Balance of Surface Strain and Reduced Dissipation 

In contrast to the momentum balance with constant /a, where the surface-strain tensor 
B simply plays no role, we now show that in the energy balance it does play a role but is 
always self-balanced. This interesting point has never been noticed before. 

The conventional energy balance reads 

P T>t (§M 2 ) =pf- u + #P + V (T •«)-$, (7) 

where 

$ = 2 faDijiTij + pSij ) = Ad 2 +2 uDijDij (8a) 

is the dissipation function. Using an identity of Truesdell, 10 we found 

$ = (A + 2p)6 2 + puj 2 - V • (2/xB • u). (86) 

Therefore, substituting T = T - 2pB into (7), the terms containing B are precisely 
canceled. This result implies: 

Theorem. For Newtonian fluid with constant shear viscosity, the work rate done by the 
viscous resistance to surface strain is locally balanced by its own dissipation rate. 

Like the skin-friction, the viscous resistance of a fluid-surface element to its strain 
usually does negative work, and hence its dissipation, —V • (2/zB -u) = § s , say. is positive. 
But, this may not be always so. 

Now. due to the theorem as well as (6), for constant pi, the balance of both momentum 
and energy of a Newtonian fluid can always be replaced by that of a fictitious fluid which 
has no surface strain at all. Indeed, the energy balance now reads 

Pjft (9 l“| 2 ) ~ Pf ‘ u + dp - V • (itll + pw x ti) - $, (9a) 
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or, for any material volume V , 



where 

$ = (A + 2 -|- hu? 2 > 0 (10) 

is the reduced dissipation , and 


t = n T = — Iln + fiw xn 

is the reduced stress due to compressing and shearing processes only. These equations can 
also be directly derived from the inner product of (6b) and tz. 

Care must be taken when providing interpretations for Eqs. (9)-( 10). We find that 
there is an exact cancelation between the portion of the dissipation and that of the work 
rate done by surface stress. Although computation efforts can be reduced by using the 
reduced stress, it would be misleading to take the “reduced dissipation” as the “total 
dissipation”. Otherwise, for example, even a solid- like rotation (say, occurred near the 
axis of a vortex) would have a dissipation. This is of course not the real physics. Like the 
case of stress tensor, the true dissipation for incompressible fluid is still given by 

$ = $ + with $ = = 2/zV • (tz • Vtz). 

With the above explanation, our results (9) and (10) imply another significant sim- 
plification and some further physical insight. While it is known that for unbounded flow 
or under periodic boundary conditions the total dissipation in V can be written as 

f $dV = f [(A + 2/i)# 2 + (iu 2 ]dV = [ (11) 

Jv Jv Jv 

so that the integration volume can be reduced to the regions with nonzero u> and d. For 
example, for a Rankine vortex, the total dissipation can be computed based on (7a), 
which requires integration of 2 fiDijD,j over the whole irrotational domain, but, by (10) 
it becomes an integration of the constant ji/w 2 over the solid vortex core. Now we have 
shown that, with constant // and the reduced dissipation function, the boundary-condition 
dependence of (11) can be removed. 

Based on the fact that for incompressible flow the viscous force depends on the 
skew- symmetric part of Vu but the dissipation rate (8a) depends on its symmetric part, 
Largestrom 13 noticed the existence of potential flow with viscous dissipation. He illus- 
trated the situation by a steady axisymmetric flow driven by a rotating cylinder of radius 
a and tangent velocity V (equivalent to the solid core of a Rankine vortex). The flow 
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is irrotational but with a local dissipation rate $ = 4 p.V 2 a 2 /r 4 and a total dissipation 
4n pV 2 . We now further see that this dissipation is exclusively from the surface strain, 
and is solely balance by the work rate done by the resistance to that strain, which is in 
turn provided by the rotating cylinder with a total power input Atv/jV 2 = f&o 2 ■ ~ ra 2 to 
keep the flow steady. Note that on a rotating solid boundary with angular velocity W . 
like the above cylinder surface, the stress t s due to surface strain can always be absorbed 
into the shear stress r. as long as the vorticity u> in (5b) is replaced by a relative vorticity 
u T = u) — 2 W.‘ Therefore, the work down by the cylinder surface to retain the flow can be 
explained as from either the surface-stain stress t s or the shear stress r. 

It should be emphasized that the stress t 9 due to surface strain may appear should we 
go beyond the balance of momentum and energy. 6,14 Its most important appearance is in 
the stress balance on an open fluid surface element. Specifically, on a fluid-fluid interface or 
free surface t 9 plays a crucial role: 8 its normal component joins the balance with pressure 
and surface tension, while its tangent components are solely responsible for balancing the 
interfacial tangent vorticity. Note, however, this appearance of t s is relevant only to the 
dynamic boundary condition on such an interface only; thus , Eqs. (6) and (9a,b) still hold. 
As a result, saving is still achieved even with an interface as the flow boundary. 

As noted already, eqs (6) and (9a,b) do not apply to flow with variable shear viscosity, 
where a direct coupling between momentum balance and thermodynamics adds additional 
terms to the S-H tensor potential (6), which, although can always be obtained, do not have 
the above neat form. 

4. Turbulent Stress Tensor versus Turbulent Force 

We now turn to turbulent flow. We first show how to achieve a computational saving 
for a given model of turbulent stress tensor, say t;j = —ufUTj = Tj{. Here and below, unless 
stated differently, u is referred to fluctuating velocity, and so are its derivatives. The 
overline means ensemble average. Once again, in the mean turbulent momentum equation 
and energy equation [or filtered equation if large-eddy simulations (LES) is considered], 
what counts is only the turbulent force r l] t = fj. say; so we may similarly replace t xj bv a 
three-component S-H tensor potential fj j, which amounts to finding the S-H decomposition 
of /. 

As an illustration, consider an unbounded or spatially periodic flow so that in the 
Fourier space spanned by the wave vector k (where we use ct,/?, 7,... = 1,2,3 to denote 
Cartesian components), there is 


f gik) — ik a T a p(k), 1 — \/ 1, 


( 12 ) 


and (1) becomes 


f 3 (k) = i[kg<t>(k) + e 3Q1 k Q A^(k)\. 


(13) 



To obtain <p(k) and A-y(k) from the given T a p(k), note that the vector identity in physical 
space 


V 2 / = V(V • /) - v X ( V X /) 


implies 


k fp(k) — kpk a f a (k) d" €f)at'Y£y6pkatkpf(i(k^. 


Comparing this and (13), it follows that 

im= k°Uk)' lA ^ = (ySp tsMtl, 

Therefore, from (12) we obtain the reduced stress tensor 

T a /3{k) = k /? “h ko t k-yT-yfiik') kpk-y Tya ( ^)] • 

It can be easily checked that 


(14) 


tkc'Totpi'k') zk q-t Q $(A?) — y)j(Ar). 


By using (14), any models of the turbulent stress tensor r aj g(Jb), linear 15,16 or nonlinear 17-21 , 
can be reduced to a S-H tensor potential before being put into numerical computation. As 
noted earlier, in principle this procedure can be applied to any bounded domain as long 
as the HWD basis therein has been established. 

Although (14) can already bring great computational saving, however, it does not 
simplify theoretical analysis since one still has to model the full Reynolds stress T a p. We 
thus propose a more thorough approach, which would lead to a significant theoretical 
simplification in turbulence modeling as well. The basic idea is that as long as an S- 
H decomposition of / is obtained, it is sufficient to directly model the three-component 
turbulent force, rather than the much more complicated six-component turbulent stress 
tensor. 

We assume the flow is incompressible for simplicity. The turbulent force can be 
expressed in terms of the “vorticity form”: 22 


i 

t> 

i 

ii 

1 

i> 

i 

ii 

(15) 

A = ^M 2 and / = u x ti 

(16a, b) 


are the mean turbulent kinetic energy and mean turbulent Lamb vector, respectively. 
Obviously, I\ is a part of the scalar potential of /, and the other part comes from the 
curl-free longitudinal part of 1. Therefore, as long as we have split 1 into a longitudinal 
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part, / = Vy, and a transverse part, l±, then we immediately arrive at the desired 

decomposition of turbulent force: 


/= -V(K + x )-l±, (17) 

which sharpens (15). Moreover, since for incompressible flow there is 

V 1 = -V 2 h^ = -V 2 A\ 

where ho is the stagnation enthalpy, it follows that if we define 

ij> = X + K, (13) 

then 

V 2 ^ = 0. (19) 

of which the solution depends on boundary values of \ and K only. Once these values are 
given, the internal values of K in the flow domain can be inferred from that of Therefore, 
the problem of modeling turbulent force exclusively amounts to modeling the mean turbulent 
Lamb vector. This observation further confirms that the Lamb vector is the key in nonlinear 
fluid dynamics including anisotropic turbulence; a rational mathematical reduction is often 
associated with a sharper physical insight. 

5. Reynolds Stress Transport versus Lamb Vector Transport 

The second-order closure models (full Reynolds stress transport) represents the highest 
level of closure currently feasible in practical Reynolds-average computations. In prin- 
ciple, these models account for more turbulence physics than lower- level models, e.g., 
two-equation models. As reviewed by Speziale 16 , full Reynolds-stress transport may be 
used as a starting point to deduce improved two-equations models under the equilibrium 
limit of homogeneous turbulence. 

The previous section clearly demonstrated that the problem of Reynolds-stress mod- 
eling can be cast to the modeling of turbulent Lamb vector. Thus, the first step towards 
directly modeling turbulent force at the level of second-order closure is to derive the trans- 
port equation for /, which is much simpler than that of the full Reynolds stress. 

From the exact Navier-Stokes equation and vorticity equation (the notations now 
stand for instantaneous quantities), it follows that the Lamb-vector transport equation 
reads 

dl 

-+u-Vl + l-Yu-uY 2 l = Q, (20 a) 


where 


Q = Yh 0 x u — 2 vu^ x t*,j 


( 206 ) 



is the source. The first term of (20b) is an inviscid coupling between longitudinal and 
transverse parts, but the second does not represent the dissipation of Z at all. We found 
that 


-2vuj x it,/ = -2i/V • (a? x D) + ^(p + = ( V x u) x w, 


\- 1 ) 


where $ = i^J 2 is the reduced dissipation. 

In order to gain more understanding of the evolution of Lamb vector, let us consider 
its square, characterized by J = l 2 j 2. It is easily shown that the J-equation reads 

d J 

+ ti • V J + Z D • Z — i/V 2 J = — i/VZ : VZ + Q j , (22a) 

where : means twice contraction so that vVl : V/ is the J-dissipation, and 

Qj = Z • Q = {'Hu) — u> 2 u) • V/io — — (Vu? 2 ■ Vu 2 + 4u> • D • D • w — 4u> • D • VTi ), (226) 

where Ti. = U3 ■ u is the helicity. In (22b), all viscous terms are from the second term of 
(20b). For two-dimensional flow Q j reduces to 

Qj = -ca 2 u • V/i 0 - ^Vw 2 • Vw 2 . (23) 

Note that for potential flow or Beltrami flow both (20a, b) and (22a, b) become the trivial 
identity 0 = 0. 

It is now clear that the structure of (20) and (22) is precisely the same as that of 
the two-dimensional equations for vorticity gradient and its square , but subjected to the 
forcing terms. Therefore, the existing studies on the latter (e.g., Novikov 23 ) can well be 
utilized to explore the behavior of (20) and (22). In particular, the increase of J is due to 
the shrinking of the fluid element rather than stretching. If the stretching is dominating 
in a flow field, the integrated “Lamb-enstrophy ,, J must exponentially decrease by this 
mechanism. This is consistent with the fact that as the turbulent eddies becomes smaller 
and smaller due to stretching, the fluctuating Lamb vector is also reduced — eventually 
we have Z ~ 0 for fine-scale turbulence, which then becomes approximately homogeneous 
and isotropic. Note that in (22b), the term -2vuj l DijDjkOJk is associated with vortex 
stretching, which also provides a negative source (a sink) of J in a highly stretched vortical 
flow, thus further reduces the Lamb vector in fine-scale turbulence. Consequently, in the 
turbulent force there remains the direct balance between the turbulent kinetic energy K 
and reduced dissipation ? = fiui 2 , which leads to the Kolmogorov law. 

To investigate the effects of the mean flow on turbulence, we perform a Reynolds 
average. The mean-flow part is now denoted by capital letters and fluctuating part is still 
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presented by lower-case letters. After subtracting mean-flow terms, the desired transport 
equation for the mean fluctuating Lamb vector reads 

31 

— + U ■ V/ + 1 ■ VU - vV 2 l = 9i + 02 + q 3 (24) 

On the right of (24). the first source is 


q i =Vho xl} — 2vLJ i x u (25) 

which has the same structure as (20b) and involves second moments at the same point. 

The second source reads 

q 2 = — (u • VZ t + V • Vu), (26) 

which contains the interaction between the mean flow and some second moments. Here, 

L + = w x U + SI x u 

is a “quasi-Lamb vector” in the sense that its first term is the analogy of Coriolis force 
acted on the mean flow caused by fluctuating vorticity, and the second term is that acted 
on fluctuating flow caused by the mean vorticity. Finally, the third source 


q 3 = -(u • V/ + / • Vu) (27) 

is a triple correlation and surely needs modeling. Equation (24) indicates that 1 is advected 
and stretched-turned by the mean flow field. 

The structure of q± clearly implies that the pressure- strain term in common modeling 
is now reduced to the correlation ofVh-o and u. Because the variation of ho is expected to 
be milder than h , and since only the regions with u> ^ 0 needs be modeled, it is expected 
that some troublesome aspects of pressure-strain modeling 16 could be bypassed. The 
second term of (25), having been identified not as the dissipation, needs some attention. 
Although it is a viscous effect, it might not be completely negligible since the derivatives 
of fluctuating vorticity and velocity are involved. A further analysis of q x could be made 
by looking at (22). 

On the other hand, we see in (24) that q 2 , the triple correlation, and the advecting- 
stretching-turning of l by the mean flow, have the same theoretical structure since they 
come from the single root; but their roles are very different. A further study of this type 
of structure is desirable. 

Obviously, various second moments are not equally important, and one can only select 
the most relevant one to derive its transport equation and close the turbulence modeling 
thereon. This is the basic spirit of second-order closure. In the usual approach, one takes 
the transport equation of Reynolds stress, in which some other second moments, such as 
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the pressure-stain correlation and dissipation-rate correlation, have to be modeled. Much 
efforts have been paid to deal with some redundant quantities and complexity because the 
full Reynolds stress is reducible. Now the tensorial Reynolds-stress equation is reduced to 
the vectorial mean Lamb-vector equation, in which the remaining quantities can no longer 
be further reduced. In this sense, we may call the above approach as an irreducible second- 
order closure. An obvious feature of this irreducibility is that the terms to be modeled are 
nonzero only in those regions where w^O; or, approximately, the regions with high peak 
of fluctuating vorticity. This is precisely the essence of turbulence as randomly stretched 
vortices. 


6. Conclusion 

In this theoretical paper we systematically explored the concept and application of 
reduced stress tensor and dissipation function, showing that they may bring significant 
savings in analysis and computation. The development is solely based on the classic Stokes- 
Helmholtz decomposition, or its modern sharpening, the helical- wave decomposition. 

In particular, for turbulent flows we propose that the study and modeling of Reynolds 
stress tensor could be reduced to that of the mean turbulent Lamb vector alone, of which 
the governing equation is closely similar to that for the two-dimensional vorticity gradient. 
This equation would be a new basis of irreducible second-order closure schemes. A numer- 
ical examination of the budget and spectra of terms in this equation, including those that 
need be modeled, is being undertaken and will be reported separately. 

References 

1. Moses, H.E., “Eigenfunctions of the curl operator, rotationally invariant He lm holtz 

theorem, and applications to electromagnetic theory and fluid mechanics,” 
SIAM J. Appl. Math. 21, 114 (1971). 

2. Lesieur, M., Turbulence in Fluids, (Kluwer, Dordrecht, 1990). 

3. Melander, M.V. and Hussain, F., “Polarized vorticity dynamics on a vortex column,” 

Phys. Fluids A 5, 1992-2003, 1993; Virk, D., Melander, M.V., and Hussain, 
F., “Dynamics of a polarized vortex ring,” J. Fluid Mech. 260, 23-55, 1994. 

4. Yoshida, Z. and Giga, Y., “Remarks on spectra of operator rot,” Math. Zeit. 204, 

235-245, 1990. 

5. Wu, J.Z., Wu, X.H., Wu, D.D., and Zhou, Y., “Intrinsic decomposition of vector field 

in bounded domain with application to fluid dynamics,” to be submitted as 
an ICASE report. 

6. Wu, J.Z. and Wu. J.M., “New aspects of the Stokes-Helmholtz decomposition and 

application to fluid dynamics,” UTSI Rept. 92-04, 1992. 


- 10 - 



7. Wu, J.Z. and Wu, “Interactions between a solid surface and a viscous compressible 

flow field,” J. Fluid Mech. 254, 183 (1993). 

8. Wu, J.Z., “A theory of three-dimensional interfacial vorticitv dynamics,” Phys. Fluids 

7. 1995, 2375-2395. 

9. Dishington. R.H., “Rate of surface-strain tensor,” Am. J. Phys. 33, 827 (1965). 

10. Truesdell, C., Kinematics of Vorticity , Indiana Univ. Press, 1954. 

11. Batchelor, G.K., An Introduction to Fluid Dynamics. Cambridge University Press, 

1967. 

12. Eraslan, R. et ah, ORNL/TM-8401, 1984. 

13. Lagerstrom. P.A., “Laminar flow theory,” Theory of Laminar Flows, ed. F.K. Moore, 

Princeton Univ. Press, 1964, 20-285. 

14. War, J.Z. and Wu, J.M., “Vorticity dynamics on boundaries,” Adv. Appl. Mech. 32, 

eds. J.W. Hutchinson and T.Y. Wu, Academic Press, 1996, 113-269. 

15. Smagorinsky, J., “General circulation experiments with the primitive equations,” Mon. 

Weather Rev., 93. 99 (1963). 

16. Speziale, C.G., “A review of Reynolds stress models for turbulent shear flows,” ICASE 

Report 95-15 (1995). 

17. Yoshizawa, A., “Statistical analysis of the deviation of the Reynolds stress from its 

eddy viscosity representation,” Phys. Fluids , 27, 1377, (1984). 

18. Speziale, C.G., “On nonlinear K — l and K — e models of turbulence,” J. Fluid Mech., 

178, 459 (1987). 

19. Rubinstein. R. and Barton, J.M., “Nonlinear Reynolds stress model and the renormal- 

ization group,” Phys. Fluids A, 2, 1472 (1990). 

20. Zhou, Y.. Vahala, G., and Thangam, S., “Development of a recursive RNG based 

turbulence model.” Phys. Rev. E , 49, 5195 (1994). 

21. Lund, T.S. and Novikov, E.A., “Parameterization of subgrid-scale stress by the velocity 

gradient tensor”, Annual Research Briefs. Center for Turbulence Research, 
Stanford Univ. (1992). 

22. Lele, S.K., “Vorticity form of turbulence transport equations,” Phys. Fluids A4, 1767 

(1992). 

23. Novikov, E.A., “Statistical balance of the vorticity gradient in two-dimensional turbu- 

lence,” Phys. Lett. A 162, 385 (1992). 


- 11 - 







REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway. Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

March 1996 Contractor Report 


4. TITLE AND SUBTITLE 

REDUCED STRESS TENSOR AND DISSIPATION 
AND THE TRANSPORT OF LAMB VECTOR 

5. FUNDING NUMBERS 

C NAS1-19480 
WU 505-90-52-01 

6. AUTHOR(S) 
Jie-Zhi Wu 
Ye Zhou 
Jian-Ming Wu 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Institute for Computer Applications in Science and Engineering 
Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23681-0001 

#. PERFORMING ORGANIZATION 
REPORT NUMBER 

ICASE Report No. 96-21 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681-0001 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA CR-198307 
ICASE Report No. 96-21 

11. SUPPLEMENTARY NOTES 


Langley Technical Monitor: Dennis M. Bushnell 
Final Report 

Submitted to Physics of Fluids. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

12b. DISTRIBUTION CODE 

Unclassified-Unlimited 


Subject Category 34 



13. ABSTRACT (Maximum 200 words) 

We develop a methodology to ensure that the stress tensor, regardless of its number of independent components, 
can be reduced to an exactly equivalent one which has the same number of independent components as the surface 
force. It is applicable to the momentum balance if the shear viscosity is constant. A direct application of this 
method to the energy balance also leads to a reduction of the dissipation rate of kinetic energy. Following this 
procedure, significant saving in analysis and computation may be achieved. For turbulent flows, this strategy 
immediately implies that a given Reynolds stress model can always be replaced by a reduced one before putting it 
into computation. Furthermore, we show how the modeling of Reynolds stress tensor can be reduced to that of the 
mean turbulent Lamb vector alone, which is much simpler. As a first step of this alternative modeling development, 
we derive the governing equations for the Lamb vector and its square. These equations form a basis of new second- 
order closure schemes and. we believe, should be favorably compared to that of traditional Reynolds stress transport 
equation. 


14. SUBJECT TERMS 

Turbulence Modeling; Reduced Stress Tensor and Dissipation; Lamb Vector Transport 


15. NUMBER OF PAGES 

13 


16. PRICE CODE 

A03 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


16. SECURITY CLASSIFICATION! 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION! 
OF ABSTRACT 


20. LIMITATION 
OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298(Rev. 2-89) 

Prescribed by ANSI Std Z39-18 
298-102 






















